Tema 5 - Flux estacionari de fluids en conductes tancats 


6 de febrer de 2020 


Introducció: 


Donat una canonada de secció circular a un moment donat coneixem sa pressió en el punt 1 i en el punt 2. 


h 


z 


Tenim que l’eix 
horitzontal 
és l’eix z 


Condició 
Re < 2300 



(Laminar) 


Si el fluid circula d’esquerra a dreta, vol dir que P\ > P 2 , la pressió a la secció 1 és major que la secció 2. 
Recordam que: 


P = 


F 

8 


Com hipòtesis direm que la pressió és la mateixa en tota la secció: 
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Que no depèn d’on estigui col·locat el punt de la pressió. 

S’objectiu del tema 5 sobretot, l’ho fonamental és trobar el perfil de velocitats. 

Que és el perfil de velocitats? 

v=o 

V=V ma\ 

v=o 

Si dibuixam un eix de coordenades on l’horitzontal és Z i la vertical és l’eix de coordenades r, veim que 
quan r és màxim, la velocitat és zero (r = màx => V 7 = 0. aplicant condició de no lliscament), i quan r és 
zero la velocitat (1/) és màxima. 

La longitud dels vectors de velocitat en funció de r han d’ésser simètrics. 


V—Vmax 


Trobar això implica resoldre les equacions del tema 3, no podem plantejar aquest experiment amb el 
procediment del tema 4, perquè hauríem de fer experiments i posar molts de sensors per saber quina velocitat 
és mou cada una de les capes del fluid. 

El tema 4 esta pensat per problemes a nivell macro, però no micro. 

Aleshores hem d’agafar les equacions del tema 1, 2, i 3 per resoldre aquest problema. 

En aquest cas: 

Canonada de secció circular 
Velocitat —> V =? 


v=o 



v=o 



Sa qüestió “és senzilla”. 

Si ho agafam amb en Newton, direm amb un exemple: seria el tir parabòlic, quina és la trajectòria d’un 
tir parabòlic? 
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Amb un sòlid. (Cas senzill amb Newton). 

A pesar d’ésser un problema senzill, sa resolució matemàtica no és senzilla (s’ha de tenir agilitat). 


1 Règims en funció del nombre de Reynolds: 

El nombre de Reynolds (Re) és un nombre que va ésser inventat per Reynolds, (Osborne Reynolds (1842- 
1912)) les seves dimensions son 1, és a dir; que no te dimensions (és adimensional). 

I ho definim com: 


Re= P^v 

í» 


p-L-V 

í» 


On: 


p : Densitat del fluid. 


D : Diàmetre del conducte. 


V : Velocitat del fluid. 


p, : Viscositat del fluid. 


L : Longitud característica del problema, en aquest cas és el diàmetre (D) 



El perfil de velocitats no depèn de la longitud (L) però si de la distància entre els punts de Vo, V max , Vq. 
Reynolds va descobrir que el seu nombre (Re) apareixia a molts d’experiments de la mecànica de fluids. 


1.1 Definicions de quan un fluid circula dins un conducte: 

Com és pot caracteritzar? 


• Fluxe és: 
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— Estacionari: Tots els paràmetres que defineixen les condicions de contorn son invariants amb els 
temps. Si no hi ha res del contorn que depengui del temps, és equivalent amb fer una fotografia a 
un temps (ti) i un altre a (to) i veim que observam el mateix. Perquè el seu contorn no depèn del 
temps. 

— No estacionari: Transitori. Tindrem un flux transitori ò no estacionari si els paràmetres si que 
varien amb el temps, els paràmetres de contorn. 


A l’exemple que havíem proposat. 



La condició que posarem és que: 


Pi = A'i P 2 = K -2 (Son constants) 

Aleshores com que les condicions de contorn son constants, el diàmetre no varia, per tant aquest exemple 
serà estacionari. 

No ha res que varii en el temps. 

Per altre banda direm que: 

• El flux és classifica: 

— Flux compressible: Els que sa poden comprimir, s’observen canvis apreciables de la seva densitat. 
p = ™ (típic a fluxos a altes velocitats). 

* Nombre de Mach: 

Si M > 0.3 és compressible. 

M = - > 0.3 

C 

(v és la velocitat del so en el medi en que es troba i c és la velocitat del so a l’aire (~ 340 [ m /s])). 
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— Flux incompressible: No existeixen canvis apreciables de la densitat. 


Que és un flux? 


Un fluid en moviment, que a travessa una superfície determinada. 


Ara podem afegir un altre classificació, en funció del nombre de Reynolds: 
• Si Reynolds: 


- Re < 2300: 

Flux laminar. Les capes, cada una d’elles va a una certa velocitat i no es creuen entre elles. 



- 2300 < Re < 4000: 


Període de transició cap a un flux turbulent. 



— Re > 4000: 

Flux turbulent, absolutament caòtic, no sa pot predir sa forma que prenen els remolins (teoria 
del caos). 





•- 

Turbuí·nl 




Com sa defineix Reynolds per una secció circular? 


Re = 


p-D-V 
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Reynolds és un indicador per quin regim estem fent feina per un cas. 


Re < 


- Laminar 

- Transitori 

- Turbulent 


En rexperiment que hem plantejat des de el principi d’aquest tema, suposàvem que Reynolds és menor 
de 2300, per tant feim feina en regim laminar. 


2 Flux laminar en conducte circular. 


Quin és l’objectiu? 

Conèixer el perfil de velocitats d’aquest flux, que ja sabem que serà laminar. 

Re < 2300 


Aquesta solució sempre serà valida per aquestes condicions sinó no serà valid. 
Per què es mou el fluid d’esquerra a dreta? 

Com a conseqüència de ses forces de pressions. 

Pl>P2 


Suposam que el flux és estacionari i incompressible. 



Quina és aquesta velocitat màxima? 
Com la podem trobar? 


V ? 

v m.n.x • 


Però el que més ens interessa, és trobar el vector velocitat () : 
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(Pe) Pj /; Pz) UjVjW 


També observam que aquest fluid tindrà una viscositat, per tant haurem de tenir en compte que és donen 
esforços tangencials. 

Primera equació per fluids per Newton. 

Recordam que l’esforç a damunt la cara y en direcció x és: 


T~yx M 


dV x 

dy 


Tema 1 


Per una banda haurem d’aplicar aquesta equació com a definició de fluid del tema 1 per saber el tensor 
d’esforços que en algun moment haurem de posar. 


y 

u 



Àm — H 


-H 


- *V 

iC 

v- —- 


v - >W 

i— 


^V—Vmax 


v=o 



V=Vmax 


v=o 


Perfil de velocitats laminar (Re < 2300). 


Si ara ens fitxam l’esforç a damunt la cara r en direcció z, podem agafar la formula anterior i aplicar un 
canvi de coordenades, per què? 

Perquè ara tenim coordenades cilíndriques. 

T v* = P • íf 

4 


r = u- 

1 rz H> dr 


Començam a aplicar equacions: 


T rz M 


dV z 

dr 


Per altre banda tenim: 


dp 

dt 


+ • {^p • = 0 (Tema 3) 


Per què? 

Perquè és l’equació de continuïtat. 
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És conservació de la massa. 

En tot problema de Newton s’ha de conservar la massa. 

En tot problema de fluids s’ha de conservar la massa. 

És una equació aplicable sempre. Sempre que no hagi reaccions químiques. 
També recordam, que si el fluid era incompressible i estacionari: 


^ + 

dt 


0 




= 0 




= 0 


p pot sortir defora del parèntesis perquè és una constant. 


x-(Ah) 


= 0 


Equació de continuïtat —> ^ ^)=°H Per estacionari i incompressible | 


És una altre equació que hem de tenir en ment: 

(V. h) = o (Equació convergent) 


Perquè la massa s’ha de conservar. 

Sa coneix com sa divergència del camp de velocitats que és zero. 

Divergència =v • Vector (Estrictament matemàtic) 


^V = l·v x + l·v v + l·v í = o 

= dx ■ V x + dy ■ V y + dz ■ V z = 0 


Així si és en coordenades cartesianes. 

Però resoldre aquest exemple en coordenades cartesianes és molt complexa. 

Aleshores farem us d’un altre sistema de coordenades. Donat que aquest sistema presenta simetria cilín¬ 
drica, ens introduïm dins les coordenades cilíndriques. 

Son polars si no hi ha z. 

Conversió de coordenades: 




La conversió la feim mitjançant trigonometria: 
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r = \Jx 2 + y 2 


7 

x = r ■ cos (ep) 

Cartesiana =>• Polar < 

Lp = arctan (J) 

Polar =>• Cartesiana < 

y = r ■ sin (</?) 


z = z 

\_ 


z = z 

\_ 


Això que vol dir? 

Que qualsevol vector a l’espai és pot expressar com: 



Cartesia Polar 

S * V ^ 

És la mateixa informació 
Qualsevol camp de velocitats és pot escriure d’aquesta manera: 

(y r , V<p, V z ) Que vol dir la velocitat amb cada una d’aquestes components? 



V z =$> El fluid és mou d’esquerra a dreta. 


V r =$> El fluid és mou d’amunt a baix. 
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K, => El fluid esta fent voltes el voltant de l’eix 


Teníem per equació de continuïtat. 


f•?)=0 


El que ens interessa és escriure aquesta equació però en sistema de coordenades polars (cilíndriques). 
Nosaltres sabem que l’operador Nabla és: 

+ + fc'jen cartesianes. 


dx 


Però hem de treballar en coordenades cilíndriques. 


En operador Nabla =v Això no é 


és correcta. 


Realment l’operador Nabla en coordenades cilíndriques és: 




..... — . i . r . JL i . JL JL 

cilíndriques y r 1 Q r ’ r dtp 5 dz 


Aleshores 


^T=(í·r.|:,i. i |4).(l/ r ,V.K) = 0 

= - • AEd-P + f . _|_ flAkl = 0 (En coordenades esfèriques) 

r or r dtp dz ^ / 


Fent un petit incís l’operador Nabla aplicat a un vector en coordenades cilíndriques 


^ _ I . O(r-JV) , 1 . 0(F V ) , Q(F X ) 

r dr ' r dtp ' dz 


Continuació: Si és vera que s’ha de complir això, igualam a zero. 
Quines conclusions podem arribar? 

El vector velocitat te components: 


^ = iy r ,v^v z ) 


Però donades les condicions de contorn: 


P\ — ki, P‘2 — 
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No varia en tota la superfície, el mateix l’hi passa a P 2 


Pl- 


I que el fluid va d’esquerra a dreta. 

Hi ha algun motiu per pensar que el fluid esta donant voltes? 



V ? 

Vtp. 

No, no hi ha cap motiu. 

Quin motiu hi haurà si hagués gir? 

A la vegada que és fa pressió, 



s’agafa un disc i és fa rotar. Aleshores si que V v tindrà un valor. Però no és el cas. 

Per tant com que no hi ha cap condició de contorn que faci girar el fluid, per condició de contorn (C.C.). 

V^ = 0 

f = (v r ,y£v)j 

I una altre conclusió que podem treure a partir de les condicion de contorn és: Sa velocitat amb la direcció 
z no depèn si estem a una posició z\ o una posició z 2 . 
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^ serà diferent de zero però no depèn de la posició 

Recordam que cada una d’elles pot dependre de les altres coordenades. 


K ( r , tp, z ) 

/ \io depèn de z 

Vz [r,(p,jn 
V<p (r, (p, z) 

El temps no apareix, perquè hem dit que era estacionari. 





K (r, ip, z) 





V z (r, (p, z) 


\ 




/ 



Anulam perquè les pressions estan fitxades i en mig no hi ha. 

= 0 


Zi Z 2 

X - X 


Una condició de contorn que l’hi faci anar més ràpid o més lent, la seva derivada parcial respecta a z a 
d’ésser igual a zero, no depèn d’ella mateixa. 


Ü§1 = 0 


Aleshores de la equació queda: 




x° 

1 _ d(,r-V r ) , 1 _ djVjtf . dWf = q 
r dr ' r ymp ' /òz 
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£ÍM = 0 

ü(f 


djVz) _ 


> Per condicions de contorn. No hi ha res que faci pensar que el fluid estigui girant ni vagi més 


dz 


= 0 


ràpid o més lent. 

Per tant ens queda: 


^ . ï? = i . ça + 1 . o + o = 

r or r 


^ ■ V = ± • = o 

r or 


Com que sabem que serà diferent a zero ^ 0) i que el producte ha d’ésser igual a zero, per tant: 


9(r-Vr) = Q 
dr 


1 

r 


d (r • V r ) 
dr 

" -V-- 

0 


On r és un paràmetre. 


Si 




això implica que r • V r és una constant. 


(r ■V r ) = K 


Aleshores=^ V r = — = K ■ - 

' r r 


Això és físicament possible o no. 

Que vol dir això? Que la velocitat de la component r (que és aquella que puja o baixe) és inversament 
proporcional a r. 

K = I< ■ - 

Que passaria si V r (r = 0) (quan r = 0) ? 

V r (r = 0) = K • ^ = oo (Físicament impossible) 

Trobam un infinit amb una condició de contorn. 

No és possible, aleshores no tenim més ramell que dir que V r és igual a una constant, però sa constant ha 
d’ésser zero perquè tengui una solució realista. 

K = 0 
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V r = 0 • 1 

' r 


Resumint: 


Ens queda: 




V r = 0 

v v = o 


|> .Justificam que els dos son zero. 


V z r, 


o \o 


Paràmetres (nomes sobreviu r) 



ip = 0 —> No gira 
V z No depèn de la posició de z 


V z (r) 


V z = V z (r) 


Que és el mateix que 



(0, 0, V z (r)) 


Si sa massa s’ha de conservar, el fluid no pot pujar ni baixar. 



Coordenades cilíndriques 


Camp de velocitats 

o 

ll 

tt> 

Per condicions de contorn 

^ = (0, 0, V z (r)) 




Per tant nomes sobreviu una component. 

Fins ara nomes hem aplicat l’equació de continuïtat. Però també en el tema 3 varem veure la conservació 
del moment lineal en forma diferencial que estava relacionat amb les equacions de Navier-Stokes. 


Tema 3 

fí = -g + P·g + ^·n 

> ^ 


Moment lineal de forma diferencial 

<§ = .p + p.t + ^ 






En aquest exercici igual que a Newton amb el tir parabòlic, hauríem d’imposar que s’ha de conservar sa 
massa i sa conservar el moment lineal. 

També sabem que: 


{F = m-lt] [Ç = p] 

Si ara ajuntam aquestes dues equacions: 

* € = -^·P + í>· 9 > + ^·rJ' 

* àÈ = dm-~à _ _ d\^ 

dv dv l·' dt 


€ = ~^ .p + p .f + ^ 
. dP _ dm-lï _ _ dP 

dv dv V dt 


Substituïm 0- 

dv 

P·W = -^ -P + p-t + ^ -T. 

Substituïm 0- 

dv 

p ~=-^.p+ P .^+f. 


No és nomes una equació, son 3, una per cada direcció. 
Idò això comença a ser una cosa més complexa. 

De fet comença a aparèixer el tensor d’esforços. 


T 

T 

T 

T 

T 

T 

T 

T 

T 


On ara aquí posarem el subíndex de les cilíndriques. 


T 

T 

r) 

Direcció r 

T 

T 

r 

Direcció tp 

T 

T 

T ) 

—ï Direcció z 

í r 

1 rr 

^~ípr 

T zr ^ 

— y T r 

T~r<p 

T W 

7~Z(p 

Tçp 

Trz 

Tpz 

Tzz j 

' T z 


Si ens fitxam ara amb el terme de l’esquerra: 
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En a que és igual el terme de l’esquerra. 


I acceleració a que és igual? 


P ■ — sabem que és a l’acceleració. 


p.dt 

P dt 


(acceleració total) 


(^ T) T 


Temal 

Es la derivada total i esta relacionada amb la part local i amb la part convectiva. 

ll (Total) = terme local + terme convectiu 

La part local (4^-) no és més que derivar la velocitat respecte el temps de manera parcial, però com que 
és estacionari és zero. (les derivades temporals son zero). 


€ = » 


La part convectiva veim que és: 




(0,0,14) 


1 d 1 d d 

r dr’ r dtp' dz 


Ja ho hem justificat O. Nabla en coordenades cil. 

(f T)T= [o- i. 1 + 0 - i·l + K-I] -(0,0,14) 

-tyf = [0 + 0 + U·A].(0,0,V’J 

(^T)T = V»·|:·U 


(o.o, v; 
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Però sabem que la variació de la velocitat (14) en relació a 2 és zero (condicions de contorn). 
Per tant: 




= v z • — 1 

2 dz 


(^ ■ 7 )' 7 = V, - 0 
(^ .?).? = 0 


Per tant el terme de l’esquerra de la conservació del moment lineal és: 


Zero igual en el terme de la dreta: 


# = -^· J P + p·t + ^·T 

• 0 = — ^·P + p·~^’ + '^·Tj 

0 = — ^ • P + p • -Tj 


Anam a veure si el terme de la dreta és igual a zero. 


= — ^ • P + p ■ ^ ' 


Descompondrem cada un dels seus termes. 

Començam amb el primer terme: 

• P =>- Recordem que estem en coordenades cilíndriques, hem d’escriure l’operador Nabla en cilíndriques 
a damunt una funció (P). 

Nabla a damunt una funció (genèric) cilíndriques: 


^ • P = ^ • 




On < (p > son vectors unitaris de cada eix de coordenades. 
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Aquest terme d’aquí novament te tres components. Una en direcció r, en direcció <p i un altre en direcció 


Per condició de contorn: 

-Sa pressió si que varia si ens movem en direcció 

z. Per tant ^ ■'z ^ 0. No la descartam. 

Pl 

P2 

X - 

— X 

[lOOOPa] 

[2Pa] 


-Però sa pressió P en cap moment deim que varia en funció de les altres coordenades. 
^•f=0No varia amb el radi (més amunt o més a baix). 

• (p — 0 No varia d’acord amb l’angle. 


Aleshores nomes sobreviu una contribució: 

^ • P = (PJr, + + %-z 

0 0 


Segon terme: 




p _ dP 

' r ~ dz 


Novament per condicions de contorn. 


P‘~Í 


4 


p ■ (9x,g y ,9z) 

g x ,g z Sa gravetat no actua ni en direcció x ni en direcció 2 . 


9y 


No juga cap paper. 


p- (0,0,0) 
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El que volem dir és que, el fluid que esta circulant per aquest conducte, no sa veu sotmès a sa gravetat. 


Perquè és mou en direcció 2 (és ortogonal a la direcció z). 
En direcció r no hi ha cap moviment . 

9-J 


I ara anam en el tercer terme: 




L’operador Nabla per el tensor d’esforços en coordenades cilíndriques. 
Aquí no tenim més ramell que escriure tot. 








\ 


T r r 

7~(fr 

T zr 


—> r 



T~Z(f 


<P 

T rz 

T cpz 

T Z z 

J 

-)■ z 


r zz => El fluid s’esta traccionant. 



Aquí novament per condicions de contorn nomes sobreviu r r 
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/ 




T, 


tpr 


T 7 , 


\ 


7~rtp ^~tptp T Z(fi 


T r 


T l 


tpz 


T z 


= < 


í ^ • T r = 
^ ■T lf = 
^ ■ T z = 


1 _ Or 
r dr 

1 Or 
r dr 


_1 . _ 0 _ 

r dtp 

_1 . _ 0 _ 

r dtp 


_ 0 _ 

dz 

_o 

dz 


1 Or _i_l _ _0_i 0_ 

r dr ' r dtp ' dz 


\T rr + 7~tpr 
\j~rtp H~~ 7~tptp 




= 


[r rz + T(p Z + r^J — 


1 dr i 1 d i_0_ 

r ' dr ' Trr ' r ' dtp ' T( ^ r dz ' T ^ r 

1 dr i 1 _0_ _i_0_ _ 

r * dr TrL P ' r ' dtp ' T W ' Oz ’ 


i 1 Or i 1 _0_ i 0 

- ’ Or ' Tr2: r ’ 0<p * T( ^ Oz * 



• r = 

■ £ = 
z = 


r Or 
1 d(rr r(f ) 


+ r 


+ 


Per tant: 


+ è 


1 _ dQyz) 

r 0(/? 


dr 

d{rT rz ) 

dr 


dtp 

1 Q( r y?y?) 

dtp 


_|_ 1 . a(ry, z ) 


0(/p 


djTzz) 

dz 


Z = 


_ 1 d(rr rz ) 


dr 


dz 

dpzy,) 

dz 

_ d(r zz ) 
dz 


■r — if 


2 : = 


^ • r z = 



1 # djrTrr) 
r dr 


+ è- 


O(Tyr) 

dtp 


+ 


1 4 O(r·Try) | 
r dr r dtp 


1 Q(r yy ) 


d(r zr ) 

dz 

d{r Z(p ) 


i 1 # ^(r yg ) 

r 0(/? 


dz 

dtjzz) 

dz 


T rz p 0 Nomes depèn de la viscositat. 

En resum que ens queda de la equació original: 


r = lf 

• £ = lf 


' _ 1 d(rr rz ) 

r dr 


0 = — ^ ■ P + p ■ ^ • Tj 


( 


= ^·Pt = ^· 


T, 


0 0 0 
0 0 0 
0 0 


\ 


0 = 


dP 

dz 


■z + o + -è • 


d(rr rz ) 

dr 


OP C i 1 0(rr rz ) 
dz * ^ * r ' Or 


Z 


Això és el que s’ha de resoldre: 
Passam a l’esquerra. 


OP 1 d(rr rz ) _ ^ 

Oz r Or 


És a dir: Derivar la pressió respecta a z és igual en el terme de la dreta (£ • ° (r £ z> )· 
Això que vol dir? 

Això vol dir que les forces de pressió sa compensen amb les forces viscoses. 
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Pi haurà d’ésser més elevada com més viscos és el fluid (haurà d’ésser més elevada sa diferència). 

I el seu moviment dependrà de la diferència de pressions que hi ha (gradient (ff ))• 

dependrà de /i (viscositat). 

Matemàticament aquí ens adonam compte que derivar una funció respecta z (terme de l’esquerra) i després 
derivam respecta a r (terme de la dreta). 

8P_ _ 1 _ d(XTrz) 
dz r dr 

Clar com que son derivades parcials perquè això sigui valid i sigui físicament possible les dues han d’ésser 
iguals a una constant (K). 

dP_ _ 1 _ d(rr rz ) _ 


dz 


dr 


Per què? 

Perquè les dues bandes derivam respecta a coordenades diferents. 
I això nomes pot tenir una solució possible si: 


dP 

dz 


K i 


1 . d(jTrz) _ 
r dr 


Idò direm que tant el terme de l’esquerra com el terme de la dreta son constants i l’hi direm K\. 
Anem a resoldre sa part de la dreta: 


i 

r 


d{rr rz ) 

dr 


= Ki 


1. Passam a la dreta r multiplicant. 


djrTrz) 

dr 


= Ki - r 


2. Ara integram tant a l’esquerra com a la dreta: (Integram respecta a r). 

■ dr = j Ki ■ r ■ dr 


L’integral respecta a una derivada sa cancel·len. 
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rr rz = f Ki ■ r ■ dr 
rr rz = Ki ■ f r ■ dr 
rr rz = K 1 ·Ç + K 2 

Hem d’afegir una constant d’integració (K 2 ): 

VTrz = Kl • g + K 2 

Ara per trobar sa constant d’integració, aplicarem condicions de contorn: 
r = 0 

• Si r = 0: 


TT rz = K 1 ·Ç+K 2 


7~rz 


K r ^+K 2 


r 


T = Ki ■ — + ^ 
1 rz IY 1 2 ' r 

t = Ki • - + Aa 

‘rz - f '-l 2^0 

A-z = ^frAl+ qf 


r r 


- Aa. 


r r ^ = oo <^= És impossible!! 


Hauria d’ésser oo però com que no és possible aleshores r rz = 0 si r = 0. 


Idò aquesta constant (K 2 = 0). Per condicions de contorn. 


Aleshores aquesta constant d’integració que hem hagut de posar, no la tenim i nomes sobreviu: 


r rz = I< l ·l 
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Que vol dir? 

L’esforç T rz és igual a una constant per r dividit entre dos: 


/ \ 7 y r dVz 

r rz (r) = K 1 · j = fi- — 

Temal 

I ara ens fitxam amb una definició del tema 1; Qualsevol esforç per Newton és pot posar proporcional a 
sa viscositat. 


fi • 


dr 


t 


Si és un fluid 
Newtonia 

I ara que una altre equació diferencial a resoldre: 

Com que per l’equació de continuïtat nomes ens hem quedat amb un terme i aquest terme ja el tenim a 
una equació com incògnita, ja nomes ho hem de justificar-ho: 

Tenim que: 


r rz = Ki ■ | 
i 


r = u -^ 

Irz P dr 


> K x • | = fi • 


dr 


K x ■ \ = fi 


dr 


Passam la viscositat a l’esquerra. 


Si ara giram l’expressió. 


r 



dV z 

dr 




dr 


Passam el diferencial de r a l’esquerra multiplicant. 


dr • ATi • ^ = h!4 

Reordenam l’expressió: 
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Ara aplicam integral. 


dV z = ïf 1 ■ r ■ dr 

z 2-fl 


I dV * = I fe ' r ' dr 

v * = fe -fr-dr 


Afegim la constant d’integració. 


u = g • >- 2 + k 3 

Arribats aquí aplicam condicions de contorn per saber que val K 3 . 
• V z (r = R) = 0 



Si teníem que: 


V z = • r 2 + K 3 

z 4 -fl 1 ° 


Si V z {r = R) = 0 


y z = j x ·r 2 

z A-fi 


K» 


0 = ^ 

4-M 


A 2 +A 3 


-A' 3 


= t x ·R 2 

4-íi 
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K a = • R- 


És el procediment de sempre; feim una integral, deixam una constant d’integració, després aplicam 
condicions de contorn i aïllam la constant d’integració. 

Llavors K 3 en aquest cas és la constant que havíem definit com K\ => ( Ií'i = jfz)- 


Aleshores ja podem dir que és V z : 


-ïl··R 2 

4-/í 


dP 

-f^-R 2 

4-n 


T _ aC 

K, = • R 1 = —• R 1 


v, (r) = • r 2 + K 3 


V, (r) = & • r 2 + K 3 


V, M = i- ■ Ki ■ r 2 + A's 


u W = i - (§?) • r2 


■t^-R 2 

4-/X 






v * ( r ) = i' (§) • r2 


(m -r 2 

4-ji \ dz ) 


V z (r) 




4-n V dz 


4-/i V dz 


Treim factor comú: 


VM = ^·m·tr 2 -R 2 


Si ara giram l’expressió: 


l '.w = è·(-S)·( R ’- r! ) 
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És a dir; aplicam equació de continuïtat més conservació d’es moment lineal més condicions de contorn, 
en coordenades cilíndriques. 

Som capaços de saber quin és el perfil de velocitats. 

És a dir; el vector 

Hem trobat una expressió que diu el següent. 

Sa velocitat nomes te component z i ens depèn de la distància en el centre del conducte, però sempre va 
en direcció z. 



I aquesta forma, observam que te una forma parabòlica. 
A que és igual el perfil de velocitats? 


• Depèn de la viscositat (/r). 

• Depèn del valor que agafa sa pressió entre el punt 1 i el punt 2 ). 

dP P2-P1 
dz L 

• Depèn del radi del conducte. 


Aleshores si nosaltres coneixem sa viscositat (/r), el radi (r) del conducte i els valors que ens donen uns 
manòmetres i la distància (L) que hi ha entre elles, a partir d’aquests valors a nivell macro podem trobar 
quin és el perfil de velocitats. 

Per cercar la velocitat màxima: 

La velocitat màxima és quan r = 0 => V max = V (just el centre del conducte). 

Sabem que: 

UM = i'(-f)'(íí 2 -P) 

Substituïm (r = 0). 
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V W = 


1 . ( dP\ ( d 2 n2 


A-n V dz 


(.fí 2 - o 2 ) 


fM 2 


T4 (r = 0) = V max — 


4 V cte 


Aleshores: 


Si ara tenim en ment que: 


U(r) 


(-f). (j R 2 -r 2 ) 


T/ — Jt 

rmax — 4-/x * V <9z V 






f )- 2 


KM = ^- 


_qt^\ [ i r 

, aJ' 1 i ? 2 


V z (r) = K 


1 — ^ j Perfil parabòlic. 


Recordam que la velocitat és: 


^ = (0,0,14) 


Ara que tenim? 

Hem trobat a que és igual el vector 


V = (0,0,14) 


On 14 = V max • (1 - ^ 


Tal com varem dir en el tema 1 el més important a la mecànica de fluids, sempre és trobar el perfil de 
velocitats. A partir del perfil de velocitats podem trobar altres propietats: 

• Exemple: El cabal volumètric ( Q ) [*/ s ]. 


Llavors sabem que el cabal és: 


^ -dÀ 


Temal 
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n = z 

dA = 2-7 t ■ r ■ dr 


Que passa ara? 

Hem de convertir aquesta integral a coordenades cilíndriques. 
En coordenades cilíndriques el diferencial d’àrea és: 


dA = 2 • 7r • r • dr ■ 'z (dA cilíndrica) 


On 2 • 7T • r és el perímetre. 


Q = ffc 


cercle 


1 


dl 


Q = fo 


R 


1 


(dl = 2-7 t ■ r ■ dr 


2 ■ it ■ r ■ dr 


Q = 


fo 


R 


4 


Vmax ( 1 fi 2 




Vmax ' I 1 


r 

R 2 


2 • 7T • r • dr 


Vmax = Constant 
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Si treim factor comú: 


Aleshores Q: 
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Q = V max ■ 2 • 7T • R 2 ■ [l - 


Q = V max ·7r·R 2 ·[l- §] 

Q V max " 7T • /? 

Per tant el cabal que circula és: 

Ens fitxam que tt • R 2 és la superfície de la secció de la canonada. 

• Més casos que podem calcular: Velocitat mitjana (V m ) [ m /s]. 

On estava definida com: 

t f _ 2 'Vmax ^ 

V ™ ~ _ 2 ' Vmax 

Tema 1 


•[« 


Q = \ • Vma X ' TT • R 2 


V 

y m 


1 

2 


Vrr 


Velocitat mitjana, que és? 

Sabem que sa velocitat en la secció no és constant, però de manera mitjana és \ ■ V max . la meitat de la 
velocitat màxima. 

Idò quina és V max 7 


V = - ■ V 

v m 2 v 771 


T / — Vmax 

v m — 2 


V — V 

y max — y m 
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2.1 Definicions: 


Nosaltres hem vist que una pressió sa pot escriure com: 

[Pa] P = p ■ g ■ h 

On h és la columna d’aigua. 

Varem veure que una pressió sa pot trobar com: 


P = p-g-h 


D’on ve això? 

Que si teníem un punt donat: 


P. Atm. 



- P. Total 


P. tot al = Pi + P.Atm 


P.total = p ■ g ■ h + P.Atm 

A vegades parlam de pressió manomètrica: 

En aquest cas: 

P.manomètrica = p ■ g ■ h 

Per què? 

Perquè P.manomètrica = P.Absoluta — P.Atm 

S’utilitza molt sa pressió manomètrica per què ja és sabut que tot esta en contacte amb l’atmosfera, 
aleshores anam a llevar l’offset, aquesta pressió que és comuna a tots i ens quedam amb la pressió manomètrica, 
del que realment mesuren els manòmetres. 

D’aquí sabem que les unitats de pressió son pascals [Pa], però donat que això (p • g ) és una constant, sa 
pressió sovint és dona en mestres columna d’aigua [ m.c.a.]. 

Per què? 

Perquè si donam la h ja sabem que: 

31 












p-g-h 


1000 • 9.81 • h 

p ■ g =>- És el factor de conversió 
Pa m.c.a. 

Aleshores per passar de m.c.a. a Pascals, basta multiplicar per p ■ g. 

P - 9 ■ h = [Pa] 
h = [m.c.a.] 

Donat uns metres de columna d’aigua (h), tenim uns Pascals i nomes hem de multiplicar per p - g 
(1000 • 9.81) . 

P - g - h = [Pa] 

D’aquí que 1 atmosfera = 101325 [Pa] Pascals, això si ho dividim entre p ■ g. 

^ = “ = ïÜSï = 10-328 KH 

També sa pot passar amb mil·límetres de mercuri (Hg). 


101325 _ 101325 


Ph g ‘9 13600 - 9.81 


= 0.759 [m.c.Hg] ~ 760 [mm.Hg] 


Aleshores parlar de metres columna d’aigua o parlar de Pascals és el mateix. 


[Pa] p H2 o ■ g-h =^Aigua-^ m.c.a. -s- = [m.c.a] 


[Pa] pHg - g-h =>Mercuri—^ m.c.m. —> ^ a } g = [m.c.m. 


Dit això, la diferència dels dos manòmetres: 


A P = Pi — P 2 (Pèrdua de pressió) 


AP= Pi^P2 
S’ha d’escriure com => p-g-h 
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I donat la unitat en metres columna d’aigua. 

És a dir; qualsevol pressió sa pot donar com a metre columna d’aigua (■ m.c.a .). 

A P = 

P'9-h 
A P = p ■ g ■ h 

Sa pressió nomes basta donant la informació de la h. 

AP = Pi - P 2 

A P = [Pa] - [Pa] 

AP = [PaHPa] = x \ m c a \ 

p-g L J 


3 Anàlisi del flux en un conducte circular. 


• Tenim que la velocitat màxima és: 


V = Ai . (-§P) 

Vmax 4 .^ V dz ) 

• Cabal volumètric és: 

Q 2 Énaa; VT • R 

• Diferencia de pressions és el que s’anomena; pèrdua de carrega. 

P\ — P 2 = p ■ g ■ h = AP (Pèrdua de carrega) 

Tenim un fluid, que per anar d’esquerra a dreta sa l’hi ha de sotmetre una pressió {P\) que ha d’ésser superior 
a (P 2 )> perquè si no, no és mourà. 

En que podem pensar en que esta relacionada aquesta pressió, per què és mogui? Factors que poden 
influir? 

• Diàmetre del conducte (. D ). 
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Cabal (Q). 


• Densitat del fluid (p). 

• El material del conducte, rugositat (e) = [L] = (mm) 

A més rugositat (e) (imperfeccions del conducte) més elevada haurà d’ésser la pressió perquè el fluid és mogui. 
Per tant, en el final, també dependrà del nombre de Reynolds (Re), no és el mateix flux laminar, que flux 
turbulent, dependrà del cas en concret. 


Cas 1 

Laminar P t (Costa d’empenya) 


Turbulent =4> P 4- 

Cas 2 

Laminar P 


Turbulent P t (Costa d’empenya) 


Aleshores aquesta diferencia de pressions és el que s’anomena pèrdua de carrega. S’ha perdut una pressió 
com a conseqüència del fregament. 

Si recordam Bernoulli del tema 3 hi havia un terme que era: 

g ■ h,L = Pi — P 2 = p ■ g ■ h,L = APl (Pèrdua de carrega longitudinal) 

Partint d’equacions conegudes: 

Hem vist que ens interessa trobar: 

/ \ 

Perquè ho feim aixi? 

/ dP\ _ p-g-fiL _ J r-> v l _ P 1 —P 2 _ Pèrdua de carrega 

\-pï) - ~L~ - S Per que > - —— - Longitud 

p ■ g ■ h L = P 1 - P 2 

\ / 

I això és molt senzill de mesurar, és posar un manòmetre en el principi (P±) i un altre en el final de la 
canonada (P 2 ) i és mesura la longitud lineal de la canonada i és fa la divisió. 

Aleshores tenim les següents expressions: 
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* V 

* Vmax ~ 4 p V dz) 

* Q o V m ax ' Tt ' R 


( 9P\ _ p-g-hL 

* l ¥/“ — r~ 


Q 2 Vmax ' Tí " R 




n = I. # 

^ 2 4-/z 


ap l = 

P • .9 • h L 

L 


Pi — P% > 0 

7T-R 2 


« = |·£·( l t‘)''·R í 


e = HS • P?) • * • * 2 = <? 


Cabal en funció de la pèrdua de carrega 


Si ara agafam aquesta expressió: I així com és útil trobar Q en funció de la pèrdua de carrega, a vegades 
els enginyers defineixen molt be el cabal i a partir d’aquí fan el contrari, trobar la pèrdua de carrega. 

Aïllam AP/y 

Q - i - S - m -*-# 

4 


A Pl (Q) = < ^ | 8 A L (Dades edifici) 


A Pl és la incògnita quan un enginyer fa el projecta. 

Per què ho feim així? 

Perquè no podem posar dos manòmetres. Així anam a preveure la pèrdua de carrega i quan la tinguem 
(A Pl), en el tema 6 veurem quina bomba impulsora (bomba centrífuga) hem de posar. 

Si ara convertim la següent expressió: 


A p QSjj^L 

LXrL 7 r-R 4 


Q = V m ·A 
< Q = V m • 7T • R 2 > 

R = f 


p _ ftElL. 

\r\ 
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Aquesta expressió, quan és valida? Sempre? 

A D Vm ·3‘2·fU,·L 

_ D 2 _ 

- Nomes en regim laminar (Re < 2300). 
- Secció circular. 

Podem fer una altre conversió de l’expressió: 


A P L = 


Vm-32 JJrL 
D 2 


128 -Q-n-L 
-K D A 


Recordam que les pèrdues lineals també és podem posar com metres columna d’aigua: 


hL — 


A Pl 
Ph 2 0-9 


m.c.a. 


Una vegada que hem analitzat el comportament del flux laminar (Re < 2300) en un conducte circular, 
ara anam a definir un coeficient de fricció que ens permetrà analitzar el flux per qualsevol regim i qualsevol 
secció. 

Sa pot generalitzar l’anàlisi de la pèrdua de carrega, i això qui ho va fer? 

En Darcy (Darcy-Weisbach) 


• Henrry Darcy (1803-1858) 


• Julius Weisbach (1806-1871) ; Va millorar la formula. 


En Darcy va definir un coeficient de fricció que englobarà tot (Reynolds, secció). 

Va dir que aquest coeficient de fricció nomenat com factor Darcy (/) és: 

Igual a sa pèrdua com a conseqüència del tram lineal (A Pl = P\ — A) dividit entre la densitat (p) per 
la longitud (L) dividit entre el diàmetre hidràulic (D H ) i després multiplicat per la velocitat mitjana (V m ) el 
quadrat i dividit entre dos. 







El diàmetre hidràulic és ( Dh ): 

Que és quatre vegades sa secció i dividit entre el perímetre. 

_ 4 'Secció 

H Perímetre 

I aquí si que pot trobar el diàmetre hidràulic d’una secció triangular, quadrada, rectangular perquè és una 
definició. Clar si sa secció és circular aleshores: 

Dh = Diàmetre del conducte 

Amb aquesta formula ja ens és igual quin és la forma del conducte. 

El que varen fer en Darcy, en Moody i altres, va ésser posar dins del laboratori diferents trams de canonades 
amb diferents seccions, diferents materials, diferents velocitats i diferents longituds i així trobar moltes f's. 

En el principi semblava un caos, semblava que no tenia sentit res, no hi havia manera d’entendre els punts 
dels resultats dels experiments. 

Ara be el que nosaltres podem fer és, com que coneixem sa solució de regim laminar per secció singular, 
anam a substituir aquesta pèrdua de carrega lineal que ja sabem a que és igual, la substituïm a la formula de 
/: 

• Factor de fricció en un conducte circular en regim laminar (Re < 2300): 

f = A P r , 

Vm ‘32-/r • L 

r _ 'V n< 32 ·/j-X-Dh -2 

DZ-p-X-V^ 

r _ 32·fi·D}j·2 

1 ~ D^Vnr 

f 6±ppDiL 
J D^-p-Vm 

Com que és un conducte circular (Dh = D). 

n 64-/^vS' 

J - Dk-p-Vm 


A TD _ V rn ·32·p·L 

— D 2 
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r _ 64 ■/! 

J ~ D-p-Vm 

En el final trobam una cosa curiosa, que = Re (Nombre de Reynolds). 

Aleshores f = . Que és la formula de Poiseuille. 

J Re 

Formula de Poiseuille 

r _ 64 

_ J ~ Re _ 

/ té dimensió? No, perquè Reynolds no en té. 

I aquesta equació f = sa la coneix com equació de Poiseuille. 

Ara ja tenim el factor de Darcy per un cas concret. 

Si ara agafam la formula general, podem treure la pèrdua de carrega longitudinal: 



La incògnita real per un enginyer és A Pl. 



Si ens fitxam la velocitat mitjana esta relacionada amb el cabal (demanda de l’edifici) però també estirà 
en funció del diàmetre que l’enginyer trii (no és el mateix una canonada de |” que una de 4”). L’enginyer 
haurà de valorar segons el diàmetre que posi, ses pèrdues que té i si serà assumible o no. 

Que vol dir? 

Que si W m son molt grosses, seran molt cares i no serà eficient. 

Sempre s’ha de minimitzar A Pl per així posar maquines petites. 

Però que implicarà això? Diàmetres grossos. 

Si volem A P L 4=^ D H 4 


• A Pl (/, Re) Per conducte circular en regim laminar: 


A P L 


f P-L-V* 

J D h -2 



38 








D = Dh Conducte circular 


A P L 


64 _ pJYkà 
Re Dfí-2 


A P L = 


64 

Re 


D-2 


Que hem trobat? 

Hem trobat una / = |?. 

Sempre i quan Re < 2300 i la secció és circula. Si no no. 

• Factor de fricció en regim turbulent (Re > 4000): 

Però que passa si 2300 < Re < 4000 “en cara no ho sap ningú que passa”. Si que sa sap, però és molt 
complexa (es tracta de transició). 

Regim turbulent (Re > 4000): Casos limit: 

— Cas 1 =^> Conducte de parets llises, és a dir que no depèn de la rugositat (e), aleshores / nomes 
és funció del nombre de Reynolds. 

/ (Re) e ~ 0 

— Cas 2 =>- Si tenim una turbulència molt grossa (Re 1"T1"), °1 que és va observar és que el 
nombre de Reynolds és tan gran que ens és igual si és 1.000.000 que 100.000.000 (Ja nedam). 
Aleshores en aquest cas / nomes depèn de la rugositat (e). 

I concretament s’ha inventat el terme de la rugositat relativa ( e r ) per a dimensionalitzar sa ru¬ 
gositat. En el fons ens és igual els mil·límetres de la rugositat si no sa relació amb el diàmetre 
(Dh). 

Clar si e = 1 [cm] (és moltíssim) però la canonada és de D H = 2 [km] , no afectarà per res (s r JJ,). 

£ — — 

D 

Aleshores si Re ttt (Sa pot menysprear Re). 

HÇr) 

— Cas 3 = > Entre els dos casos limit, / dependrà tant per el nombre de Reynolds (Re) com per la 
rugositat relativa (e r ). 
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/ (Re, e r ) 


En el fons s’ha observat que depèn del nombre de Reynolds (Re) i de la rugositat relativa (e r ). 

Això amb els experiments que es varen fer. Però quan estem en els casos límits: 

{ - Per canonades llises. 

- Per Re ttt • 

Després van arribar una serie de gent: 

Aquí nomes n’anomenam quatre caso. 

Perquè amb aquests quatre casos més el cas de || (aquests cinc casos), així més o menys ho tenim tot. 

• Factor de fricció segons el règim. 


— Cas de Blausius: el que va trobar és que: 


f = 0.3164 • Re~ 0 25 


e r ~ 0 Canonades llises. 
3000 < Re < 100.000 


f(Re) 


— Cas de Von Karman i Prandtl: 


7F = - 2 ' 0 • l°Sio 


2.51 

Re-Vf 


s r ~ 0 Canonades llises. 
Re > 100.000 


f(Re,f ) 


Per cercar / ho hem de fer mitjançant un mètode numèric iteratiu. 

Cas de Nikuradse: 


77 = - 2 -°; lQ gio • fe] 


Canonades molt rugoses. 
Re ttt (Enorme) 


/ 0h) 


— Cas de Colebrook-White: 


7 J = -2.0 • log 


10 


s r i 2.51 
3.7 ' Re^f 
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Canonades rugoses. 
Re > 4000 


/ (e r ,Re) 


Després de tot això, va arribar un senyor anomenat Lewis Ferry Moody (1880-1953) i va unificar tots 
els resultats dels experiments dels personatges anteriors a una gràfica anomenada diagrama de Moody. 



Pèrdua de carrega lineal. 


A P L = f- 


Dh-Z 


On / ho agafam del diagrama de Moody. 

Això és el que va fer Moody, arrels dels estudis previs ho va dibuixar tot. 
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I que sa va trobar? 

A que és igual f (Re) (f en funció de Reynolds) però també depèn de la rugositat relativa (s r ). 
Aleshores Moody va fer corbes per cada valor de rugositat relativa (ey). 



— Exemples de trobar f amb el diagrama de Moody: 
Re = 100.000 ] 

1. > =* / = 0.0251 
s r = 0.002 I 

Re = 50.000 I 

2 . > =>■ / = 0.022 

s r = 0.0004 1 

• Procés per trobar / amb el diagrama de Moody: 


1. Trobar línia del valor del nombre de Reynolds (Re). 



2. Trobar línia de la rugositat relativa (e r ) 
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0.001 

0.0008 

0.0006 

0.0004 

0.0002 


0.0001 

0.000,05 


IÓ® 


3. Trobar a on creuen Re i e r i veure quin és el valor de l’eix de l’esquerra per saber /. 



Flujo completamente turbulento 


0 000,001 


0 000,005 .. 


0 08 


0 08 


0.07 


D 08 


0 M 


0 06 


0 03 


0 06 


0 02 


0.015 


0 04 


0.01 

0.008 


0.006 

0.004 


0 03 


0.025 


0 DOS 


D.02 


0.001 

0.0008 

0.0006 


e (ft) I e (mm) 
Acero remacfiado -0.01 i 3 

Concreto -0.001 10.01 0.3-3 

Madera -0.001 0.3 

Hierro colado O.OOOtó 0.26 

Hierro galvamzado 0.0005 0.15 

Hierro forjado O.OOOÉ5 0.046 

Tubos estirados 0.000b05 0.0015 

I I 1 I Mhl I llHU I H ll! 


0.0004 

0.0002 


0.015 


0.0001 

0.000.05 


0.01 


[ 008 


0 . 000.01 


0.006 


2 345679 2 345679 

10 5 106 
Número de Reynolds Re 


Seguidament podem veure uns exemples de rugositat relativa dels materials més típics que s’empren per la 
edificació com per l’enginyeria civil: 
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Material canonada 

Rugositat absoluta (e) [mm] 

“Acero Bridado” 

0.9-9 

Acer comercial 

0.45 

Acer galvanitzat 

0.15 

“Concreto” 

co 

1 

co 

o 

“Concreto Bituminoso” 

0.25 

CCP 

0.12 

Ferro forjat 

0.06 

Ferro de fos 

0.15 

Ferro dúctil 

0.25 

Ferro galvanitzat 

0.15 

Ferro dolç asfaltat 

0.12 

GRP 

0.030 

Polietilè 

0.007 

PVC 

0.0015 
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Canonada o revestiment 

Rugositat absoluta (e) [mm] 

tubs estirats d’acer 

0.0024 

Canonades de llautó o coure 

0.0015 

Fundició revestiment de ciment 

0.0024 

“Fundición con revestimiento Bituminoso” 

0.0024 

“Fundición centrifugada” 

0.003 

Fundició asfaltada 

0.06-0.18 

Fundició 

0.12-0.6 

Acer comercial soldat 

0.03 - 0.09 

Ferro forjat 

0.03 - 0.09 

Ferro galvanitzat 

0.06 - 0.24 

Fusta 

0.18-0.9 

Formigó 

0.3-3 

“Acero roblonado” 

0.9-9 


• Pèrdua de carrega lineal; 

És una pèrdua d’un tram de canonada (APl). On la pèrdua és: 


A P L = f- 


2 


/ 

Ho trobam per formules o per diagrama de Moody. 

V m 

Valors tabulats en reglament (C.T.E.), normativa (0.5 — 2 [ m /s]). 

D h 

? (No el sabem) 


Però que passa? 

Que les instal·lacions no nomes estan fetes de canonades, també estan fetes de elements singulars. 
Aleshores quina pèrdua de carrega provoquen aquests elements singulars? 

És el que s’anomena com a pèrdua de carrega localitzada. 
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4 Pèrdues de carrega localitzades (pèrdues singulars). 



Com és mesura aquesta pèrdua? 

Clar si el càlcul que hem fet per una secció lineal, ara imaginam que ficam la geometria de l’element 
singular i intentam resoldre les equacions, no ens en sortirem perquè és bastant més complexe. En el final en 
que s’ha tradueix? 

Que el fabricant esta obligat a dir quina és la pèrdua de carrega d’aquell element. 

Quina és la pèrdua de pressió que l’aigua experimenta quan circula per aquest element singular, com a 
conseqüència de la fricció, i no nomes dependrà de la rugositat si no que també dependrà de la geometria. 


Pi AP P2 


r 


n 

It-cw. ilux 


te. flux 


Cal destacar que aquestes pèrdues poden ésser molt majors a les pèrdues pròpies de les canonades (pel 
motiu que sigui una instal·lació té més colzes que trams rectes). 

Si la xarxa de la canonada és molt gran, aquestes pèrdues no és poden modelitzar de forma individual. 
Aleshores el que és fa és suma totes les pèrdues singulars que és una equivalència a una pèrdua d’una 
longitud concreta. 

Per tant es suma aquestes pèrdues a les pèrdues longitudinals. 
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Realment estem afegint longitud a l’instal·lació. 

Pèrdues longitudinals + Pèrdues singulars 

4 

Longitud total = Longitud instal·lació + Longitud lineal equivalent 

4 

Pèrdua total de l’instal·lació = Longitud instal·lació + Longitud lineal equivalent 
Anam a veure un cas concret, tenim un colze: 

Aquí tenim l’equació de Bernoulli del tema 3. 

^ + ~2 + 9 ' z i ~ w m = ^ + -2-+9' z 2 + g'hL + g·hs 

Aquesta equació el que ens diu és que si tenim un punt 1, un tram de canonada lineal fins un punt 2, la 
informació que tenim en el punt 1 i en el punt 2 és: 

Pressió (Pi, P 2 ), densitat {p\, fn). Velocitat (Ifi, V 2 ), i la cota (z \, z-x) respecta a un punt. 

Si tenim aquesta informació i ames sabem que per en mig del recorregut hi existeixen màquines, llavors el 
que ens relaciona el punt lamb el punt 2 per conservació de l’energia és aquesta equació (Bernoulli). 

El fluid que esta en el punt 1 té una energia i l’ha de conservar quan arriba en el punt 2. 



Quan sa perd energia? 

• A la interacció amb les màquines. 

(bomba ( I energia) o turbina (- energia)) 

• Pèrdues de fregament, pèrdues lineals. 

= p- g-h L 

On g-hL son les pèrdues longitudinals. 
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Però ara és que adames podem trobar elements singulars. 

Aleshores el que farem serà escriure l’equació de Bernoulli i afegir la pèrdua d’elements singulars de tal 
manera que l’energia en el punt 1 més l’energia transferida per la màquina, ha d’ésser igual a l’energia en el 
punt 2 més les pèrdues produïdes pel recorregut. 

Com és pot modelar el terme de pèrdues singulars? 

Per un cas en concret: Per un element en concret, s’ha demostrat que és pot escriure de la següent manera: 

V 2 

(Pèrdues singulars) g • hs = K • -^ L 

On K (Re, geometria) i K és un factor que obligatòriament el fabricant l’ha de donar. 

El que diu aquesta expressió és que la pèrdua de carrega de l’element singular és proporcional a l’energia 
cinètica multiplicat per un factor K. 

K és constant? 

No sempre, depèn de Reynolds i de la seva geometria. 

Aquesta expressió també és pot escriure com: 


g·h s = K· v -f 


Q = A • v m 

y = Q 
Vm A 


A = 7r • r 2 


r = — 

' 2 


D = D 


H 


g ■ h s = K ■ \ 


Q 

*_Rjl 


g • h s = K 


8 Q 2 


Una altre manera de expressar les pèrdues és mitjançant el terme de longitud equivalent (L e ). 



Quina és la longitud equivalent? 

Hi ha fabricants que donen la K o la longitud equivalent ( L e ), les dues coses son el mateix. 
Sabent que la A Pl (/) és: 
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APr = f . P±Xrn 
L J D h -2 


AP l = p- g-h L 


'R·9·hL = f·'*pt L = L e 


9 -h L = / • 


L.-Vrl 
Dh *2 


* 


També sabem que: 


9 • h s = K • Y * 


Si ara igualam els dos tipus de pèrdues: 


9 • h L = f 


g·h s = K· y -f 


Dh-2 

/2 
m 

2 


Si ara aïllam la longitud equivalent (L e ). 


9 ■ h L = g ■ h s 
f LeU = K Yà. 

J D h -2 2 


f = 

J d h \ 


K ■ 


f-k = K 



Si ara agafam g ■ = g ■ hs- 


9-hs = f 


Dh-2 


Sabem que| = Q = Q 


, aleshores: 


g· h S = f· ~D^ 2 ~ 



, , M-2,) 

g·h 3 = f· -fe 1 


g-h s = f 


*'(%) 


Dh-2 
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Per què escrivim això? 

Perquè les pèrdues singulars (g ■ hs ) son proporcionals a velocitat el quadrat i proporcionals també a cabal 


el quadrat. 

Per què? 

Perquè V m és proporcional a /i. 
Segons fabricants: 


K 


K(V m ) 
~^K(Q) 


En el final que ens queda de Bernoulli? 


Pl l V 1 l P'2 l V 2 

i t _L _ L- _L n . 'V- — m — ± z _L _ è- 


P 1 


+ 9 ' z i — w m — ^ + ^+ 9'Z2 + g-hL J rg-hs 


g ■ hL 

On { ^ son els termes afegits. 

9 • h s 


Si tenim un conjunt de elements singulars, hem de fer un sumatori de pèrdues singulars: 


9 • h s = EKi • f 


Recordam que: 


[Pa] A P s = p ■ g • h s 

Si ho volem saber amb metres columna d’aigua ( m.c.a .). 

= h s [m.c.a] 


hL = Pèrdues per fregament de l’aigua dins de la canonada/conducte. 
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